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Abstract 

In previous papers we have considered mutual simulation of n-partite pair- interaction 
Hamiltonians. We have focussed on the running time overhead of general simulations, 
while considering the required number of time steps only for special cases (decoupling 
and time-reversal). These two complexity measures differ significantly. Here we derive 
lower bounds on the number of time steps for general simulations. In particular, the 
simulation of interaction graphs with irrational spectrum requires at least n steps. We 
discuss as examples graphs that correspond to graph codes and nearest neighbor inter- 
actions in 1- and 2-dimensional lattices. In the latter case the lower bounds are almost 
tight. 

1 Introduction 

Simulating Hamiltonian evolutions of arbitrary quantum systems on a quantum com- 
puter is an idea that goes back already to Feynman [|J. This would be an attractive 
application of future quantum computers since there are no known efficient classical 
algorithms for simulating generic dynamics of many particle systems. Here a quantum 
computer can be any quantum system provided that its time evolution can be con- 
trolled in a universal way. In particular, the problem of controllability includes the 
question which Hamiltonian evolutions can be simulated efficiently by the considered 
Hamiltonian system [2-13] Assume the natural Hamiltonian H to act on an re-fold 
tensor product Hilbert space 

H n := n <8 H (g> • • • (g> H 

where each TL denotes the Hilbert space of a qudit, i.e. a d-dimensional system. The 
Lie algebra su(d) of traceless Hermitian operators on TL is a m := (d 2 — l)-dimensional 
real vector space. Let B := {a a \ a = 1, . . . , m} be an orthogonal basis of su(d) with 
respect to the trace inner product (A\B) := tr(A^B)/d for A,B £ su(d). 
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Then the most general qudit-qudit interaction on n coupled qudits is given by 



k<l a0 

where J is chosen to be a real symmetric mn x mn-matrix with zeros for k = I. Note 
that the symmetry of the coupling matrix J does not imply any physical symmetry of 
the interaction. It is a consequence of our redundant representation that turns out to 
be very useful. The coupling matrix J consists of m x m-blocks. The m x m-matrix 
Jki given by the block at position (k, I) describes the coupling between the qudits k 
and I. We have = J^, i. e. the matrix describing the coupling between the qudits I 
and k is just the transpose of the matrix describing the coupling between k and I. The 
blocks on the diagonal are zero matrices. 

In the setting discussed here and in most other articles on simulation of Hamiltoni- 
ans the only possibilities of external control are given by local unitaries on each qudit. 
We assume that one is able to implement them independently. Formally, all control 
operations are elements of the group C := SU(d) <8> SU(d) ® • • • <g> SU(d). A common 
approximation is to assume that all operations in C can be implemented arbitrarily 
fast ("fast control limit"). The simulation of Hamiltonians is based on the following 
"average Hamiltonian" approach that has successfully been used for describing Nuclear 



Magnetic Resonance experiments since many years (e.g. [14, [L5[). 

Let t\, t2, ■ ■ ■ ,t n be positive numbers and u\, 112, ■ ■ ■ , un £ C be control operations. 
Then the algorithm 

perform u\, wait t\, perform u\, perform u%, wait ti, perform u\, . . . , per- 
form un, wait t]\[, perform njy 

implements approximatively the unitary evolution 

N 

3=1 

If the times tj are small compared to the time scale of the natural evolution according 
to the natural Hamiltonian H this is approximatively the evolution according to the 
average Hamiltonian 

J2t jUj H U ]/T, (1) 
3 

where r := ^ s ^ ne s l° w down factor of teh evolution, i.e., the time overhead of 

the simulation. For investigating the time overhead and the number of time steps of 
a simulation it turns out to be useful to work with the coupling matrices instead of 
considering the Hamiltonians themselves. To express the effect of the control operations 
on the coupling matrix J note that any unitary operation u S SU(d) corresponds to a 
rotation on the m-dimensional sphere via the relation 

u ^ ( Cg^a) u = ^ c a a a , 

a a 

where the vector c = (ci, £2, . . . , c m ) is obtained by applying a rotation U £ SO(m) on 
the vector c = (ci,C2, . . . ,c m ). It is straightforward to verify that conjugation of Hj 
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by v := (g) uy> (8) ■ ■ ■ ® u( n ) corresponds to conjugation of J by a block diagonal 
matrix of the form 



v -.= u w © u {2) © • • • e u {n) e so{ 



m) 

fc=l 



The condition for correct simulation is hence given by 

j=J2w v ^' ( 2 ) 

J 

where the orthogonal matrix Vj- corresponds to the unitary Vj for j = 1, . . . , N. 

The question of optimal simulation has been completely solved so far only for the 
case of two-qubit Hamiltonians [0, 16]. Optimal simulation protocols are constructed 



in [ 16 1 . The number of time steps is between 1 and 3. 



For n qudits we have shown in M, 17, 11] that the eigenvalues of the coupling matrix 



J and J provide lower bounds on the simulation time overhead. Graph theoretical 
notions can also provide upper bounds. In certain cases the bounds are known to be 

tight § Hi- 



Here we do not focus on the time overhead but on the number of time steps in 
the n partite case. Upper bounds for decoupling (switching off the Hamiltonian) and 
time-reversal (simulating of —H by H) are directly given by the parameters of known 
schemes ]|, ||, || [Xlfl . Lower bounds for time-reversal have been derived in [17|. These 



results shown that time overhead and the number of time steps are not connected 
in any obvious way. Here we address the general problem of mutual simulation of 
pair-interactions on n qudits. 

In Section 2 we derive lower bounds on the number of time steps for the general 
problem of mutual simulation of pair- interactions. The type of coupling is assumed 
to be the same between all nodes, but the strengths and the signs may vary. Like 
the bounds on the time overhead || 0, |ll), the bound derived here make us of the 
spectrum of the corresponding coupling matrices. In the special case that one wants to 
cancel some interactions while keeping others ( "the simulation of a certain interaction 
graph") the bounds refer directly to the spectra of the corresponding adjacency matri- 
ces. Surprisingly, the greatest lower bound can be given if the smallest eigenvalue of 
the adjacency matrix is irrational. In Section 3 the bounds shall be applied to three 
different cases where long range interactions are present between all spins. The first 
two cases are the simulation of nearest neighborhood coupling in a square lattice and 
a cyclic chain. The third case is a certain interaction graph that has been proposed to 



prepare the states of a graph code [18, 19, 2C]. 



2 Lower bounds on the number of time steps 

We restrict our attention to interactions with an additional symmetry, namely Hamil- 
tonians of the following form 

k<l a/3 
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The matrix W := (wki) is a real symmetric n x n-matrix with zeros on the diagonal. 
It describes the coupling strengths and the signs of the interactions between all qudits. 
The matrix C = (c a p) is a real symmetric m x m-matrix characterizing the type of the 
coupling. This means that all qudits interact with each other via the same interaction 
and that only the coupling strengths and the signs vary. It is important that in this 
special case the coupling matrix J can be expressed as a tensor product of W and C, 
i.e., J = W®C. 

To derive a general lower bound for simulating arbitrary interactions J by a tensor 
product interaction J = W (8> C it is useful to observe that the decisive condition in 
eq. (0) is invariant with respect to the following rescaling of interactions: Multiply 
each m x m-block k, I of J and J with the same factor r^. In the case W (8> C we 
can therefore assume w.l.o.g. that W is a matrix with only 1 as non-diagonal entries 
as long as we restrict our attention to Hamiltonians with complete interaction graphs, 
i.e., Wki 7^ for k ^ I. Formally, rescaling is denoted as follows. Let A/B be the 
entry-wise quotient of the matrices A and B provided that A has only zero entries at 
those positions where A has also a zero. Then we consider the simulation of J/(W®I) 
by K ® C where K is the matrix with only one as non-diagonal entries and in the 
diagonal and / is the matrix with only 1 as entries. 

In order to emphasize that the assumption Wki = 1 is not an assumption on the 
real physical coupling strength we formulate the following theorem for general W and 
use rescaling only in the proof. 

Theorem 1 (Lower bound) Let J := W ® C be the coupling matrix of the system 
Hamiltonian, J an arbitrary coupling matrix of the interaction that is simulated, and 
H the time overhead. Denote the minimal and maximal eigenvalues of C by A m ; n and 
Amax; respectively and its rank by r(C). Let I be the m x m-matrix whose all entries 
are 1. Let s be the number of eigenvalues of J/(W (g> I) that are not contained in the 
interval 

I • — [ /^max) M^min ] • 

Then the number of time steps required to simulate Hj by Hj is at least s/r(C). 



Proof: The condition for a scheme tx, V\,t%, V2, ■ ■ ■ , ijv, Vn to be a simulation of J reads 

TV 

Y,tjVj(W ® C)Vf = J. (4) 
3=1 

Since the matrices Vj are block-diagonal we can rescale each m x m-block such that 
we obtain 

N 

tjVj{K C)Vj = J/{W 8) I) . 

3=1 



In the following we denote the rescaled coupling matrix J/(W <8> I) by J' . 

:f =1 ^(i®c)^ 

JV 



Set R := Y^f =1 tjVj(l ® C)Vj '. By adding R on both sides we obtain 



^2 tjVj ({K + 1) <S> C)vf = J' + R . (5) 
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The rank of the matrix {K + 1) is 1 since all its entries are 1. Consequently, the rank 
of the left hand side of eq. (||) is at most Nr(C). 

Set a := //A m i n and b := //A max . Note that the eigenvalues of R are contained in the 
interval [a, b]. The rank of the right hand side is at least the number of eigenvalues of 
J' outside the interval X = [—b, —a]. This is seen as follows. 

Let P and Q be the projections onto the sums of all eigenspaces of J' with eigen- 
values smaller than —b and greater than —a, respectively. Clearly, s is equal to the 
dimension of the image of P ®Q. Denote by J and R the s x s-submatrices of J 1 and 
R defined by (P © Q)J'(P © Q) and (P © Q)R{P © Q), respectively. 

Due to the choice of P and Q the spectrum of J is contained in the interval 
(— oo, — b) U (— a, oo). The spectrum of R is contained in the interval [a, b] since the 
minimal (maximal) eigenvalues of a matrix cannot decrease (increase) when projecting 
the matrix. 

We prove the theorem by showing that J + R has full rank. Set / := ^(a + 6). From 
the triangle inequality it follows that for every unit vector |^) £l s one has 

\\(j+R)\m = ii(J+/i+A-/i)i*>n 

> ii(j+/i)i*>n-ii(A-/i)i*>ii. 

The eigenvalues of the shifted operators J+/1 and R—fl are contained in the intervals 
( — oo, — i(6 — a)) U (\{b— a), oo) and [— | (b — a), |(6 — a)], respectively. This implies 
that 

||(i + /i)|*>||-||(A-/i)|*>||>o 

because the norms can be bounded by the absolute values of the eigenvalues: ||( J + 
fl)\^)\\ > ||(6 - a) | and \\(R - < \\{b - a)\. This completes the proof. □ 

The following upper bound can be easily derived from Caratheodory's theorem |^] . 

Theorem 2 (Upper bound) Every simulation that is possible can be achieved within 

n(n — 1) 2 

-m 2 + 1 

2 

time steps. 

Proof If J can be simulated by J with time overhead r then J/r is in the convex 
span of the matrices Vj JVj with notation as above. The dimension of this convex set 
is at most m 2 n(n — l)/2 since the diagonal blocks are empty and each matrix VjJVj 
is symmetric. Caratheodory's theorem states that each point in an M dimensional 
convex set can be written as a convex sum of at most M + 1 extreme points. □ 

There are interesting cases where the bound of Theorem [l] can be tightened. Assume 
J = W®C with the same matrix C as the interaction that is used for the simulation. In 
other words, only the strengths and the signs of some interactions should be changed. 
By the same rescaling trick as above we can consider the problem to simulate (W /W)® 
C by K © C where K has only 1 as non-diagonal entries. 

Theorem 3 (Lower bound) Let W®C be the coupling matrix of the natural Hamil- 
tonian and W © C the coupling that we want to simulate. Assume all non-diagonal 
entries of W to be non-zero. 



5 



1. Let C be a positive semidefinite matrix. Then the number of time steps is at least 
the number of positive eigenvalues of W/W. 

2. Let C = diag(l, 1, . . . , 1) be the m x m-identity matrix. Then the number of time 
steps is at least n — k, where k is the multiplicity of the smallest eigenvalue /J, m i n 
of W/W. 

3. Let the natural coupling be C := diag(0,0, 1), i.e., we have a z ®a z interactions 
between all spin-l/2-particles. Let the set of local control operations be restricted 
to ia x -transformations. Then one requires at least n — k time steps with k as 
in Case 2. If /x m i n is irrational then at least n steps are necessary. In any case, 
n(n — l)/2 + 1 time steps are always sufficient. 

Note that Case (3) is of special interest since it deals with simulation procedures 
that do not rely on any first order approximation. In this case all summands in eq. (||) 
commute and the unitary transformation implemented by the simulation scheme coin- 
cides exactly with the exponent of the average Hamiltonian. 

Proof (of Theorem |3|) 

Case 1 This statement is a corollary of Theorem [l] since J/(W ® I) = (W/W) <S> C. 
The number of positive eigenvalues of this tensor product matrix is r(C) times the 
number of positive eigenvalues of W/W if r(C) is the rank of C. 

Case 2 Consider the rescaled problem to simulate A ® C by K <g> C. In this case the 
right hand side of eq. (||) reduces to 

A®C + rl, (6) 

where r := ^ • tj is the time overhead and A := W/W. In jy|] we have shown that r 
is at least — fx m ™- The reason is that the spectrum of tK has to majorize the spectrum 
of A since A is a convex sum of conjugates of tK. Therefore the rank of the matrix in 
expression @ is at least m(n — k). Since the left hand side of eq. (||) has at most the 
rank mN the number of time steps is at least n — k. 

Case 3 In this case one can use a more convenient formalism. We drop the matrix C 
and characterize the interactions by W and W. If a zz-interaction is conjugated by 
iir^-transformations only two possibilities occur. The interaction term between spin k 
and spin I acquires a minus sign if exactly one spin of both is subjected to a conjugation. 
The interaction is unchanged if either no spin or both spins are conjugated. Instead of 
representing the time step j by the 3n x 3n block diagonal matrix Vj we can represent 
it by an n x n diagonal matrix Xj. The diagonal entries are ±1 and indicate which 
spins are subjected to conjugation. Then eq. (Q) reduces to 

N 

Y,hX 3 KX 3 =A. 
We add the identity matrix on both sides. Due to XjlXj = 1 we obtain 

j j 
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The rank of the left hand side is at most the number N of time steps. To estimate the 
rank of the right hand side note that the time overhead r := ■ tj is at least — /M m i n . 
Hence the dimension of the kernel of A + rl can at most be the multiplicity of the 
eigenvalue /x m j n . This shows that the number of time steps is at least n — k. 

Assume fJ, m in to be irrational. Then the time overhead r is necessarily greater than 
—fi m in- This can be seen as follows. The optimization with respect to the time overhead 
reduces to the following convex problem. Consider the matrix XjKXj for an arbitrary 
time step j. Its non-diagonal entries are ±1 and indicate which interactions acquire 
a minus sign in jth step. In graph-theoretical language, the set of matrices that can 
occur as XjKXj are exactly the Seidel matrices of complete bipartite graphs (see ||). 
Then the optimal r is the minimal positive number such that A/t is in the convex span 
of the set of Seidel matrices of complete bipartite graphs. Geometrically, the convex 
span is a polytope having the Seidel matrices as its extreme points. It is embedded in 
the n(n — l)/2 dimensional vector space of real symmetric matrices with zeros on the 
diagonal. Let O be the origin. Consider the semi-line vA for v > 0. Then the optimal 
simulation is the unique intersection point P of the semi-line with the boundary of the 
polytope. The quotient of the distance between and A and between and P is the 
optimal time overhead. This quotient can never be irrational. The reason is that P has 
rational entries since it is the solution of a linear equation over the field Q of rational 
numbers. Hence r is greater than —fi m m and A + rl has necessarily full rank. This 
proves that we need at least n time steps. 

To see that n(n — l)/2 + 1 time steps are always sufficient we can argue as in the 
proof of Theorem |l] with Caratheodory's theorem. The dimension of the convex span 
of the matrices XjKXj is at most n(n — l)/2. □ 

In the following section we will consider the task to cancel some interactions and 
keep others. Then A := W /W has only 1 and as entries. In graph-theoretical 
language, it is the adjacency matrix of the desired interaction graph. 

It is surprising that it is relevant for our lower bound whether the smallest eigenvalue 
of the adjacency matrix is irrational (see Case 3.). It is not clear whether this is only a 
feature of our proof or whether there is a true connection to the irrationality of graph 
spectra. 

In the following section we will use graph spectra for deriving lower bounds based 
on the above theorems. We show that some of them are almost tight by sketching 
simulation schemes based on well-known results on selective decoupling. 

3 Applications 

In this section we consider the simulation of special interaction graphs that appeared in 
the literature in various applications. Some interesting models in quantum information 
theory refer to quite idealized types of Hamiltonians like nearest neighbor interactions. 
If the natural Hamiltonian contains long range interactions between all nodes one may 
try to simulate the idealized interaction. Then the problem is to cancel the unwanted 
terms without destroying the desired interactions. The examples below show that it 
may cause a large number of time steps to cancel unwanted long-range interactions 
no matter how fast they are decreasing with the distance. As long as they are not 
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Figure 1: Open and cyclic spin chains 

neglected, the control sequences that cancel them may be rather long. Here we restrict 
our attention to a z (g> o z interactions between n qubits. 

3.1 Cyclic and open spin chains 

We examine first a quantum system of n spins equally spaced on a one-dimensional 
lattice. The interaction between the spins may decrease with the distance between 
the spins. We consider the problem to simulate a chain with only nearest neighbor 
interactions, i.e. we have to cancel the interactions between all non-adjacent pairs. 



The desired interaction graphs can be seen in Fig. 3.1 . 

The adjacency matrix of the cyclic spin chain is S + S T , where S denotes the cyclic 
shift in ra dimensions. The eigenvalues of S and S T are the rath roots of unity. The 
eigenvalues of S + S T are given by 

23fi(exp(i27rZ/n)) with 0</<ra. 

For odd n the lower bound on the number of time steps is n since the least eigenvalue 
is irrational for all n > 3. If n is even it is rational and has multiplicity 1. Hence the 
bound is n — 1 in this case. There are numbers ra where this bound is almost tight. 
This is shown by the following example. Let ra be an even number with the property 
that a Hadamard matrix of dimension n/2 exists. This is for instance the case for each 
power of 2 (see p2|| ). We construct a simulation scheme that consists of 2 subroutines. 

The first subroutine simulates the interaction between the pairs {1, 2}, {3, 4}, . . . , 
{ra — 1, n} and the second simulates {ra, 1, }, {2, 3}, . . . , {ra — 2, n — 1}. Note that all 
the pairs in the same subroutine are disjoint. The problem to simulate the interactions 



between disjoint pairs is a special instance of well-known "cluster decoupling" [11 1 
where the interaction between independent cliques are cancelled and the interactions 
within the same clique remains. It can be achieved using Hadamard matrices having 
at least the number of cliques as dimension (compare 0, The entries ±1 in column 
j determine which spins are conjugated in the step j. The dimension is the number of 
time steps of the decoupling subroutine. 

Using this method, we need ra/2 steps in each subroutine. Therefore we have given 
a simulation with ra steps, whereas our lower bound is ra — 1. In general the number 
of steps for simulating the circle grows only linearly in ra. This shows that the lower 
bound is quite good even for general ra. 



The adjacency matrix corresponding to the open spin chain has the eigenvalues [23 



2cos( — — : -i) , i = l,...,n. (7) 
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i) 

Figure 2: Simulation of the square lattice interaction with 4 subroutines 

The smallest eigenvalue is irrational for all n > 2. By Theorem |3| (Case 3) we conclude 
that the number of time steps is at least n. 




3.2 Square lattice 

We consider a quantum system of n = I 2 spins located on a two-dimensional square 
lattice. For simplicity assume I to be even. We want to simulate a lattice with only 
nearest neighbor interactions. 



The desired interaction graph is shown on the left of Fig. 3.2. This kind of interac- 
tion can for instance be used for preparing the initial state in the 'One- Way Quantum 
Computer' proposed in [24]. The eigenvalues of the corresponding adjacency matrix A 



are known in graph theory [23]: 



2cos (rri*) +2cos (rri J ) ' i >-? = 1 ---> L ( g ) 

We first consider the time overhead. An upper bound is given by 4 since this is the 
chromatic index of the graph ||. It is easy to see that the minimal eigenvalue of A is 
given by 

A min = 2cos(^/) +2cos(^-/) . (9) 

By Theorem ^ (Case 3) the lower bound on the number of time steps is n since the 
smallest eigenvalue is irrational. Note that this example shows that the complexity 
measures time overhead and number of time steps may differ significantly. 

An upper bound on the number of time steps can be obtained as follows. The graph 
has 2(1 — 1)1 edges. We can partition the edges into 4 sets of edges such that each set 



contains only disjoint interacting pairs. These 4 partitions are shown in Fig. 3.2. The 
simulation consists of 4 subroutines simulating one the interactions in one of the 4 
classes. For each subroutine we choose Hadamard matrices with a dimension that 
is at least the number of cliques. The numbers of cliques are l 2 /2 or l 2 /2 + I in each 
subroutine. Since there exist Hadamard matrices for every power of 2 the square lattice 
graph can always be simulated in 0(l 2 ) = 0(n) time steps. 



3.3 Graph codes 

The computational power of different n spin interactions is not well understood yet. It 
would be interesting to know which n qubit transformations can easily be implemented 
when a certain Hamiltonian is given. However, one of the few examples where the 
power of specific Hamiltonians is directly used (without using them to implement 2- 
qubit gates) is the preparation of states of graph codes proposed in pL9| . Here the 
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1 2 

Figure 3: Interaction required for preparing states of a graph code of length 5. 



codes states are obtained by the free time evolution according to a Hamiltonian with 
a z ® a z interactions. The graph representing a code is the interaction graph that can 
be used for preparing the states. We assume the natural interaction to be equal zz 
interaction between all 6 spins, and would like to simulate the interaction graph in 



Fig 3.3 



The eigenvalues of the 'wheel' in Fig. |3.3| can easily be computed by any computer 
algebra system. They are 

1 1 1 1 1 1 1 1 

1 + V6, 1-V6, --V5--, --V5--. 

The minimal eigenvalue is —1/2 — \/5/2. By Theorem [| (Case 3) the minimal num- 
ber of time steps is 6. An implementation with 12 time steps is given as follows. The 
scheme consists of 3 subroutines each consisting of 4 time steps. As above each sub- 
routine simulates the interaction between disjoint cliques and cancels the interaction 
between different cliques. Subroutine 1 has the cliques {1, 2, 6}, {4, 5}, {3}. The clique 
partitions in subroutine 2 and 3 are {3, 4, 6}, {1, 5}, {2} and {1}, {4}, {2, 3}, {5, 6}, re- 
spectively. In each subroutine decoupling the different cliques can be achieved by 
Hadamard matrices of dimension 4 since no subroutine has more than 4 cliques. Hence 
we have 4 time steps in each subroutine. 



4 Conclusions 

We have derived lower bounds on the number of time steps for simulating arbitrary 
pair-interactions between n qudits. Like the lower bounds on the time overhead, they 
make use of the spectrum of the coupling matrices. However, there is no direct con- 
nection between both complexity measures since the time overhead refers to spectral 
majorization while the bounds on the number of time steps refer to the number of 
eigenvalues not contained in a certain interval. We have shown an example where the 
number of time steps is of the order n but the time overhead is independent of n. 
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